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The  conditions  leading  to  intense  low 
frequency  generation  and  strong 
localization  in  two  level  systems 


Yuri  Dakhnovskii1  and  Horia  Metiu 
Department  of  Chemistry  2  and  Physics 
University  of  California 
Santa  Barbara,  CA  93106 


Abstract:  A  charged  particle  in  a  double  well  can  emit,  when  driven  by  an 
intense  laser,  intense  low  frequency  radiation  or  become  strongly  localized  in  one 
of  the  wells.  By  using  perturbation  theory  and  asymptotic  analysis  we  determine 
the  conditions  for  which  these  processes  take  place,  for  a  two  level  model. 
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The  behavior  of  a  charged  particle  in  a  symmetric  double  well  driven  by  a 
strong  laser  has  been  studied  recently  in  several  articles^1*7!.  Grossman,  Dittrich, 
Jung  and  Hanggit1*2!  have  shown  that  the  laser  can  prevent  a  particle  located 
initially  in  one  of  the  wells  from  tunneling  to  the  other  well.  Bavli  and  Metiul8-9J 
demonstrated  that  a  semi-infinite  Gaussian  pulse  can  take  the  electron  from  a 
delocalized  energy  eigenstate  of  the  double  well,  localize  it  in  one  of  the  wells  and 
keep  it  there. 

The  emission  properties  of  an  electron  in  a  double  well  are  also  interesting. 
Numerical  solutions  of  the  time  dependent  Schrodinger  equation^7!  show  that  the 
Fourier  transform  jx(G)  of  the  induced  dipole  p(t)  has  peaks  at  three  types  of 
frequencies:  shifted  even  harmonics  (SEH)  (Q  =  2nco  ±  A  with  n=l,2,...),  pure  odd 
harmonics  (Q  =  (2n  +  l)co,  with  n=0,l,2,...),  and  the  frequency  Q  =  A.  The  shift  A  is 
a  function  of  the  laser  frequency  co  and  its  intensity  I.  These  parameters  can  be 
chosen  to  make  A  arbitrarily  small  and  p.(Q)  has  a  very  low  frequency  peak.  The 
presence  of  this  peak  was  called1^  low  frequency  generation  (LFG).  For  certain 
values  of  {co,I},  called^  points  of  accidental  degeneracy  (AD),  A  becomes  equal  to 
zero  and  p(t)  acquires  a  static  component.  In  a  symmetric  double  well  the 
presence  of  a  static  dipole,  induced  and  maintained  by  the  laser,  is  possible  only 
if  the  electron  density  has  a  static  part  that  is  asymmetrically  distributed  between 
the  wells,  i.e.  if  electron  localization  occurs. 

The  calculations  also  show!7!  that  when  (co,I}  approaches  an  AD  point  the 
LFG  intensity  is  sometimes  very  large,  exceeding  the  intensity  of  the  fundamental 
(i.e.  the  component  having  the  frequency  D  =  co);  as  {<o,I}  reaches  an  AD  point  the 
system  has,  sometimes,  a  very  large  static  dipole.  This  means  that  the  laser 
induced  localization  is  strong. 

Finally,  it  is  also  observed  that  if  {<o,I)  is  not  at  an  AD  point  then  p(G)  has  no 
even  harmonics  (EH)  (i.e.  no  peaks  at  D  =  2nco,  n=l,2,..),  in  agreement  with  the 
selection  rules  which  state  that  such  peaks  are  forbidden  to  all  orders  in 


perturbation  theory.  As  {co,I}  approaches  an  AD  point  and  A  ->  0,  the  SEH  peaks 
(at  the  frequencies  Q  =  2nco±A)  get  closer  arid  closer  and  become  pure  even 
harmonics;  strong  even  harmonic  generation  is  observed,  in  disagreement  with 
the  selection  rules. 

Clearly  the  shift  A  is  a  central  parameter  in  describing  this 
phenomenology;  interesting  processes  take  place  when  A  is  small  or  zero.  In 
this  article  we  use  a  combination  of  perturbation  theory  and  asymptotic  analysis 
to  provide  a  simple  formula  for  A  and  to  find  the  conditions  under  which  the  LFG 
is  intense  and  the  localization  is  strong. 

To  do  this  we  assume  that  the  processes  described  above  can  be  represented 
in  the  space  spanned  by  the  two  lowest  energy  eigenstates  1 1)  and  I  2)  of  the  bare 
(i.e.  no  radiation)  Hamiltonian.  R.  BavliDOl  validated  this  assumption  by  showing 
that  the  time  dependent  wave  function  for  the  double  well  model  used  in  Refs.  6 
and  7  can  be  represented  at  all  times,  with  better  than  90%  accuracy,  by  the  linear 
combination  a ](t) 1 1)  +  a 2(t) 1 2). 

We  use  the  two  level  Hamiltonian: 

H  *  e{  1 1X1 1  - 1 2X2  !)}*{!  1><2 1-12X11  )}^l2E(t)=  eaz  -  p12E(t)  cx.  (1) 

Here  1 2  )  and  1 1  )  are  the  ground  and  first  excited  states  of  the  electron  in  the 
double  well,  respectively,  and  ox  and  a2  are  Pauli  matrices.  The  zero  of  energy  is 
halfway  between  the  energy  levels  £i  and  £3,  and  2e-e1  •  e2. 

The  induced  dipole  is: 

n(t)  =  p12  (y,t  l  cx  I  y,t),  (2) 

where  ly,t)  satisfies  the  time  dependent  Schrodinger  equation  with  the 
Hamiltonian  (1).  By  using  standard  methods  we  derive  for  ji(t)  the  equation  of 
motion: 

dp(t)/dx  =  -(g/hfo)2  {cosCe0sin(x)]  f  dt’plT)  cos[e0sin(t’)]  + 

JO 
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with  e0  =  2p.l2  Eg/fico  and  t=o)t. 

By  usingt^11^ 

exp[i  e0sin[t]=  £  Jn(e0)exp[in'c]  (4) 

where  Jm  are  Bessel  functions  of  integer  order  we  can  rewrite  Eq.(3)  as: 
dji(t)/dt  =  -(eJ0(e0 i/hco)2  J  dx,p(t’)  -  (e/hco)2  F(r,|i)  (5) 

with: 

F(t;p)  =2  £  J2n{e0)cos[2nt]  f  dtV(T)[Jo(eo)  +2  £  J2n(e0)cos[2nT]]+ 
n  =  1  JO  k  =  l 

2  £  J2n+l(®o)  sin[(2n+l>r]  f  dt’ji(x)[  2  £  J2k+i(eo> cos  C2nt’]l  (6) 

n  =  1  JO  k  =  0 

The  equation  (6)  can  be  formally  “solved”  to  give: 

\i(x)  =  cos  [A  i/co]  -  (e/hco)2  f  dt*  cos[A(t-t’  )]  F(t’;^)  (7) 

J  0 

We  identify  the  frequency: 

A  =  (eJ0(e0)/h)  (8) 

appearing  in  the  first  term  in  the  rhs  with  the  shift  A  discussed  earlier  in  this 
article.  Therefore,  the  AD  points  are  the  points  in  the  {oo.I}  plane  given  by  2n12 
Eg/hco  =rn  where  rn  are  the  zeros  of  J0. 

In  the  zeroth  order  approximation  (i.e.  if  e/hco  ->  0)  Eq.(7)  gives 
p.(t)ssji12cos[At]  and  the  Fourier  transform  p.(D)  has  a  peak  at  the  frequency  A. 
When  A  is  small  this  corresponds  to  low  frequency  generation.  In  this 
approximation  the  intensity  of  the  LFG  peak  is  very  high,  as  compared  to  other 
peaks,  since  p.12  is  the  highest  amplitude  a  Fourier  component  can  have. 

We  must  now  examine  the  first  order  correction  p1(t)  and  try  to  answer 
several  questions.  Will  the  higher  order  terms  modify  the  frequency  A?  Will  they 
change  the  amplitude  p.12  of  the  zeroth  order  term?  Under  what  conditions  is  the 
zeroth  order  term  larger  than  all  the  others? 
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The  first  order  correction  is: 

H/t)  js  -  (e/ho))2  I*  *’  cosCt-t*]  F(t';|i0).  (9) 

This  is  obtained  by  replacing  p.(t)  under  the  integral  in  Eq.  (7)  with  ji0.  The 
analysis  of  this  term  is  tedious  and  we  present  here  only  the  conclusions.  (1)  We 
find  no  indication  that  the  first  or  the  higher  order  terms  will  renormalize  the 
shift  A.  Floquet  theory1^  shows  that  in  a  two  level  system  there  could  be  only  one 
shift  frequency^.  Therefore,  for  all  parameters  for  which  the  perturbation 
expansion  is  valid  the  shift  is  given  by  Eq.  (8).  (2)  The  first  order  correction 
contains  terms  whose  time  dependence  is  given  by  cos[At].  As  a  result,  the 
corrected  dipole  has  the  form  p.(t)  =  [pi2  -  (e/hco)2  a]  cosAt,  where  a  is  a  sum  of 
Bessel  functions  squared  and  is  always  smaller  than  one.  Thus  as  long  as 

e/hto«l  (10) 

this  amplitude  correction  is  small.  (3)  |i(1)  contains  shifted  and  pure  harmonic 
terms.  The  amplitude  of  these  terms  is  of  the  form  (e/hw)2  b  where  b  is  less  that 
one;  as  long  as  Eq.  (10)  is  satisfied  the  amplitude  of  these  terms  is  smaller  than 
p.0.  (4)  The  first  order  corrections  contain  “secular”  terms  of  the  form 

t(e/ho>)2Asin(At)  which  become  infinitely  large  in  time.  This  is  a  frequent 
nuisance  in  the  time  dependent  perturbation  theory.  It  does  not  mean  that  ji(t) 
diverges,  but  only  that  as  we  increase  t  we  are  stepping  outside  the  radius  of 
convergence  of  the  expansion.  Thus  we  are  certain  that  P12  cosAt  is  the  dominant 
contribution  in  p(t)  only  if  Eq.  (10)  and 

t  (e/Tico)2  A  «  1.  (11) 

are  satisfied.  Note  that  Eq.  (11)  implies  that  the  time  when  the  secular  terms 
become  important  is  longer  and  longer  as  A  0;  there  are  no  secular  terms  if 
A  =  0. 

The  condition  (11)  implies  that  it  is  possible  that  LFG  is  intense  and  the 
localization  is  strong  only  for  a  finite  time,  even  if  (e/ho))2  <  1.  Our  numerical 
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experiments  found  no  evidence  for  this,  which  suggests  that,  when  calculated  to 
all  orders,  the  secular  terms  add  up  to  a  well  behaved  function  whose  values  are 
smaller  than  (i0. 

To  see  under  what  conditions  the  secular  terms  will  lead  to  corrections  that 
are  smaller  than  p.0  we  have  performed  an  asymptotic  analysis  of  ji(t).  We  start 
from: 

H(x)/Hi2  =  l-(cZhco)2  J  ^  Jn<eoWe(P  f  exPk  n  *il  f  dt2  exp[-i  n  t2](i(t2)/^12 

=  —  Jo  Jo 

which  is  obtained  by  integrating  Eq.  3  then  using  Eq.  4.  This  equation  can  be 
solved  by  successive  iterations.  The  first  order  term  in  the  iteration  scheme  is 
obtained  by  putting  fa(x2)=l  under  the  integral  in  the  right  hand  side.  Performing 
the  integrals  leads  to  an  expression  of  the  form 
l-(e/hco)2  [J0(e0)2('c)2/2  +  x  d  +  g  +  fit)] 

Here  d  and  g  are  constants  smaller  than  one  and  fit)  contains  sines  and  cosines  of 
ncot,  where  n  is  an  integer.  As  the  iteration  proceeds  the  term  1  -(e/ho)2  Jo^)2  x2/2 
leads  to  cos[(A/co)t].  This  cosine  will  exceed  the  next  fastest  growing  secular  term 
if  J0(e0)/2  x  »  d.  For  the  series  leading  up  to  the  cosine  to  converge,  we  must  also 
have: 

x  (e/ho)  J0(e0)  <  1  (12) 

These  two  conditions  give  the  time  interval  in  which  cos  At  is  the  largest  secular 
term  in  the  expansion: 

1/2  « cot  J0  <  l/(e/hco)  (13) 

If  Jq  is  very  small  the  asymptotic  analysis  ensures  that  (ii2  cos[At]  is  the  dominant 
term  in  p(t)  at  long  times;  moreover,  condition  (  8)  ensures  this  dominance  for 
short  times.  Together,  these  conditions  cover  the  whole  territory.  They  state  that 
if  (e/hco)  «  1  and  if  A  is  small  we  have  intense  LFG;  if  (e/hco) «  1  and  A  =  0we  have 
good  localization. 

BavliW  tested  these  results  by  solving  the  equation  of  motion  for  |i(t) 
numerically.  In  Fig.  1  we  show  (fit)  for  (e/hco)  «  1,  when  the  predictions  made  by 
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the  theory  are  expected  to  work.  We  find  that  if  the  parameters  are  such  that 
Eq.  (8)  predicts  that  A  =  0  the  calculated  ji(t)/p12  is  equal  to  a  constant  which  is 
close  to  one,  plus  a  term  oscillating  with  a  small  amplitude  and  a  high  frequency. 
|i(t)/p12  is  equal  to  one  at  all  times  only  if  the  electron  is  perfectly  localized,  at  all 
times,  in  one  of  the  wells.  The  example  shown  in  Fig.  1  comes  very  close  to  this. 
We  found  no  case  in  which  the  predictions  of  the  theory  were  in  error.  Fig.  1  also 
shows  a  case  in  which  (e/hco)  «  1  and  the  value  of  A  given  by  Eq.  (8)  is  small.  The 
time  evolution  of  p(t)  shows  a  large-amplitude-low-frequency  component  and  a 
high-frequency-low-amplitude  one;  as  predicted  by  theory  the  term  (i12  cos(At) 
dominates  the  behavior  of  p(t)  and  LFG  is  intense.  We  have  also  calculated  the 
Fourier  transform  of  p(t)  and  obtained  from  it  the  value  of  A  and  the  absolute  value 
of  the  Fourier  component  p(A).  We  varied  the  laser  frequency  and  the  laser  field 
intensity  E0  so  that  e0=2p12Eo/ho)  =6.  The  theory  predicts  that  as  long  as  (e/hco) «  1 
Eq.  8  is  valid  and  A  is  proportional  to  (e/hco).  The  plot  shows  this  to  be  true.  The 
deviations  from  this  prediction  are  not  large  even  for  (e/hco)  of  order  2.  Fig.  2  also 
shows  the  absolute  value  of  p(A)/p12  which  is  the  intensity  of  LFG.  As  long  as 
(e/hco)«l  this  is  close  to  one  which  is  the  maximum  intensity  for  this  model.  In  a 
double  quantum  well  p12  is  rather  large  leading  thus  to  large  LFG  emission.  The 
predictions  of  the  present  theory  also  work  well  for  almost  all  the  AD  points 
discovered  numerically  in  Ref.  7.  The  few  discrepancies  occur  probably  because 
the  system  studied  in  Ref.  7  is  approximately  a  two  level  system. 

Previous  work^9’10^  has  used  perturbation  theory  to  calculate  the  energy 
difference  AE  between  the  Floquet  quasi-energies.  LFG  and  SHG  were  not 
discussed,  but  the  connection  to  localization  was  made^1,2^:  localization  can  occur 
only  if  the  Floquet  quasienergies  are  equal;  in  the  present  context  this  means 
A  =  0.  However,  previous  work  used  the  high  field  limit  and  obtained  AE  under 
more  restrictive  conditions  than  those  given  here.  The  question  of  how  large  a 
fraction  of  the  electron  density  is  localized  and  how  intense  LFG  is  has  not  been 
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addressed  before.  We  find  that  these  effects  are  large  if  (e/hca)  «  1.  While  the 
field  strength  is  important  for  making  A  small,  it  is  not  a  player  in  determining 
whether  the  LFG  is  intense  or  the  localization  is  strong. 

While  we  used  quantum  wells  as  an  example,  the  two  level  analysis 
presented  here  is  more  general  and  these  processes  may  be  detectable  for  other 
systems.  For  example,  Corkum  and  his  collaborators^12^  have  observed  them 
while  exposing  ions  to  a  strong  short  pulse. 
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Figure  Captious 

Fig.l.  The  induced  dipole  p(t)  divided  by  the  transition  dipole  |x12  as  a  function 
of  time  obtained  by  solving  the  equation  of  motion  numerically.  (1)  e0=2[i12E0/hct)  = 
2.4,  which  corresponds  to  a  point  of  accidental  degeneracy  (i.e.  A  given  by  Eq.  (8)  is 
zero),  hence  to  localization;  (2)  e0=2p12Eo/hci)  =  2.7.  The  population  oscillates  in 
accordance  with  Eq.  (7).  Since  e/hco  «1,  fi(t)  is  dominated  by  the  first  term  in 
Eq.  (7). 

Fig.2.  The  dependence  of  the  shift  frequency  and  of  the  amplitude  p(A)  of  the 
low  frequency  mode  on  e/hco  for  2y.12E(yho)  =  6.0. 
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